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The longitudinal resistances of a six-terminal graphene p-n junction under a perpendicular mag- 
netic field are investigated. Because of the chirality of the Hall edge states, the longitudinal resis- 
tances on top and bottom edges of the graphene ribbon are not equal. In the presence of suitable 
disorder, the top-edge and bottom-edge resistances well show the plateau structures in the both 
unipolar and bipolar regimes and the plateau values are determined by the Landau filling factors 
only. These plateau structures are in excellent agreement with the recent experiment. For the 
unipolar junction, the resistance plateaus emerge in the absence of impurity and they are destroyed 
by strong disorder. But for the bipolar junction, the resistances are very large without the plateau 
structures in the clean junction. The disorder can strongly reduce the resistances and leads the 
formation of the resistance plateaus, due to the mixture of the Hall edge states in virtue of the 
disorder. In addition, the size effect of the junction on the resistances is studied and some extra 
resistance plateaus are found in the long graphene junction case. This is explained by the fact that 
only part of the edge states participate in the full mixing. 
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I. Introduction 

The successful fabrication of graphene, a monolayer of 
carbon atoms arranged hexagonally, 2 ~ 4 had fueled many 
experimental and theoretical works. For undoped 2-D 
graphene sheet, Fermi energy is located at the Dirac neu- 
tral point. Around the Dirac point, the graphene has a 
linear dispersion relation, which leads to quasiparticles 
obeying the massless Dirac-like equation and presents ex- 
traordinary properties. 5-7 For example, for a graphene 
under a strong perpendicular magnetic field, its Hall 
plateaus assume half-integer values, 3,4,8 h/[g{n + l/2)e 2 ), 
where g = 4 is the spin and valley degeneracy and n 
an integer. By varying the gate voltage, both carrier 
type and concentration of the graphene sheet can be 
tuned. 9 A graphene p-n junction is formed by connect- 
ing up one p-type graphene and one n-type graphene. 
Many exciting phenomena closely related to the massless 
Dirac character of carriers, 10-14 such as relativistic Klein 
tunneling 10,11 and Veselago tensing 12 , are predicted for 
graphene p-n junctions. 

Recently, the electron transport through the graphene 
p-n or p-n-p junctions was extensively investigated both 
experimentally and theoretically. 15-19 In quantum Hall 
regime, Williams et al. ls experimentally found that the 
two-terminal conductance exhibits plateaus with half- 
integer values, g(\n\ + l/2)e 2 /h, in the case of unipolar 
junctions and fractional values for bipolar junctions. At 
about the same time, the theoretical works by Abanin 
and Levitov 19 explained the appearance of the fractional 
plateaus by means of the mixture of the electron-like and 
hole-like Hall edge states in the vicinity of the junction 
boundary. There were some subsequent investigations on 
graphene junctions. Using Anderson short-range disorder 
potentials, Long et al. 20 and Li et al. 21 numerically com- 



puted and analysed the transport behavior of graphene 
junctions. They found that conductance plateaus emerge 
in the case of suitable disorder strength. Also, T. Low 22 
considered the long-range interface and edge disorders 
in the armchair, zigzag, and antizigzag edge graphene 
ribbons, and numerically simulated the result of the con- 
ductance plateaus. 

Very recently, Lohmann et al. 17 measured the Hall and 
longitudinal resistances in a six-terminal graphene junc- 
tion device. In their experiment, the difference of carrier 
concentrations between two adjacent regions (the left and 
right regions) is introduced by chemical doping, and a 
global gate voltage controls the carrier concentrations in 
the two regions. By tuning the gate voltage and dop- 
ing densities, the graphene ribbon can be of p-p, p-n, 
n-p, or n-n type. They found that Hall resistances of 
the left and right regions exhibit half-integer plateaus, as 
usual. Furthermore, the longitudinal resistances also ex- 
hibit plateau structures. In particular, the longitudinal 
resistances at opposite edges are not equal. For instance, 
with the filling factors {vl,vr) — (2,-2) (here vl/r is 
the Landau filling factor in the left /right region), the 
longitudinal resistance at one edge is h/e 2 but it is zero 
at the other edge. By using the concept of the mixture 
of Hall edge states near the p-n junction boundary, they 
explained the appearance of these plateaus and the differ- 
ence between the longitudinal resistances at the opposite 
edges. So far, there is not any theoretical investigation 
which gives quantitative and numerical result to explain 
Lohmann's experimental result. More effort needs to be 
done in order to find out how disorders affect the resis- 
tance plateaus and how the plateaus depend on disorder 
strength. 

In this paper, we theoretically and numerically study 
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electron transport through graphene junctions. Fol- 
lowing the experiment by Lohmann et al, we con- 
sider the six-terminal graphene junction device [shown 
in Fig. 1(a)]. A perpendicular magnetic field B is ap- 
plied to the graphene sheet. By using the tight-binding 
Hamiltonian and the Landaucr-Buttikcr formulism in the 
framework of non-equilibrium Green's function method, 
the longitudinal and Hall resistances are calculated. The 
numerical results show that the longitudinal resistances 
R t and Rb of the top and bottom edges, respectively, are 
usually different in the presence of the magnetic field, as 
expected from the property of the chirality of Hall edge 
states. There is an essential difference in the transport 
behavior between unipolar and bipolar graphene junc- 
tions. For unipolar (n-n and p-p) junctions, the longi- 
tudinal resistances have plateau structures in the case 
of clean (no disorder) graphene. The resistance plateaus 
can keep in the moderate disorder strength, and they are 
destroyed until the very strong disorder. On the other 
hand, for bipolar (n-p and p-n) graphene junctions, the 
longitudinal resistances Rt and Rb are very large in the 
clean case. But in the presence of disorder, they are 
strongly decreased even when the disorder strength is 
weak. For n-p junctions, the top-edge resistance Rt re- 
duces to a moderate positive value, but the bottom-edge 
resistance Rb drops to zero or even turns negative. In 
a suitable range of disorder strength, the resistances Rt 
and Rb of bipolar junctions also have plateau structures 
due to the full mixing of Hall edge states. For the low- 
est filling factors, the resistance plateaus exist in a very 
broad range of disorder strength. Hence, they are pro- 
duced easily in experiment. But for high filling factors, 
the plateau only emerges in a narrow disorder range, if 
it exists. Furthermore, for moderate disorder, the plots 
of the longitudinal resistances versus the gate voltage ex- 
hibit plateau structures in both cases of unipolar and 
bipolar junctions and the plateau values are only deter- 
mined by the filling factors (vl, vr), which are in excel- 
lent agreement with the recent experiment by Lohmann 
et al. In addition, we also find some extra resistance 
plateaus in long graphene junctions. This is explained 
by the suggestion that only part of edge states partici- 
pate in the mixing mechanism. 

The rest of the paper is organized as follows. In Section 
II, we describe the model and give the details of our cal- 
culations. In Section III, the numerical results are given. 
Finally, a brief conclusion is presented in Section IV. 

II. Model and calculation 

We consider the six-terminal graphene junction device 
which is illustrated in Fig. 1(a). The six terminals are 
labelled by the terminal-z (i = 1,2,..., 6). Terminals 1 
and 4 are the current source and drain, respectively, and 
other terminals are used as the voltage probes. The whole 
graphene device is basically divided into three regions: 
the left, central and right regions as shown in Fig. 1(a). 
Each of the left and right regions includes three terminals. 
The dimension of the central region is described by the 
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FIG. 1: (Color online) The panel (a) is the schematic dia- 
gram for the six-terminal graphene p-n junction device and 
the panel (b) is for the Dirac energy in the left, central and 
right regions. The panels (c) and (d) show the resistances vs. 
Er at the clean graphene junction (W = 0) with El = —0.1 
(c) and —0.2 (d) , respectively. The other parameters are 
the width TV = 50, the length M = 20, and cj> = 0.01. 



integers N and M [see the red-site (dark gray) region 
in Fig. 1(a)]. There are totally 2N x (2M + 1) carbon 
atoms in the central region. The width and length of the 
central region are a(3N— 1) and y/3a(M +1), respectively. 
The width of voltage terminals is chosen as V3aN : where 
a « 0.142 nm is the distance between two neighboring 
carbons. Fig. 1(a) shows the case N = 3 and M = 10. 

The whole system is subjected to a perpendicular mag- 
netic field which leads to the formation of Landau levels. 
In quantum Hall regime, bulk states are compressible, 
and the chiral edge states flow along the edges. This 
behavior is independent of the type of the edges if the 
graphene ribbon is wide enough. So we choose wide 
zigzag-graphene ribbon for our simulation. Compared to 
the hopping term, both Zeeman splitting and spin-orbit 
coupling are very small and hence they are negligible. 
Furthermore, we adopt Anderson on-site disorder. In 
fact, other types of disorders (e.g. the interface disorder 
and the long-range disorder) may exist in the experimen- 
tal device and also lead to the edge state mixing. But 
from the point of the edge state mixing, these other types 
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of disorders should have the similar effect. Here we can 
assume Anderson disorder only exists in the central re- 
gion based on the fact that the effect of disorder in the 
left and right regions is suppressed by edge states when 
the ribbon is wide enough. 23 

In the tight-binding representation, the Hamiltonian 
of the six-terminal graphene junction is given by 20,21 ' 24 



H = eia\a,i 



where a] and a% are respectively the creation and annihi- 
lation operators at site i, and tt the energy of Dirac point 
(i.e., the on-site energy). In the left and right regions, ej 
is equal to El and En, respectively [Fig. 1(b)], which can 
be controlled by the gate voltage in the experiment. In 
the central region, e, = k(E R - E L )/(2M + 2) + E L + «J», 
where the column index k = 1, 2, 2M+1 [see Fig. 1(a)] 
and Wi is the on-site disorder energy. W{ is assumed 
to be randomly distributed in the range [—W/2,W/2], 
where W is the disorder strength. The second term in 
the Hamiltonian stands for the nearest-neighbor hopping. 
The effect of the magnetic field B is addressed by the 
phase (pij = J? A-dl/<fro in the hopping interaction Hamil- 
tonian where A = (—By, 0,0) is the vector potential and 
4>o = file.. The magnetic field B is applied to the whole 
device (including the six terminals and central region) 
along the perpendicular direction. 

The multi-terminal resistance 25 is defined as Rtj t ki — 
(Vk — Vi)/Ii^j, where the contacts i and j are terminals 
used to draw and input current, and the two contacts 
k and I are used to measure the voltage difference. We 
introduce two longitudinal resistances Rt = i?i4,23 (on 
the top edge) and Rb = i?i4,65 (on the bottom edge) and 
two Hall resistances Ri = Ru } 26 ( m the left region) and 
R r = i?i4 i 35 (in the right region). These four resistances 
obey the relation, 



<ij> 



(te^alaj + h.c), (1) 



Rt + R r = Rb 4- Ri ■ 



(2) 



From the Landauer-Biittiker formula at zero temper- 
ature, the current flowing into terminal-i is given by 
h = (2e 2 A)£ j(7ti) Tij{E F )(yi - V,) 26 Here, T l3 (E F ) = 
Tr[riG r r,-G a ] = 1,2,..., 6 and i ^ j) is the trans- 
mission coefficient from terminal- j to terminal-z at Fermi 
energy E F , Ti(E F ) = i[Sl(E F ) - H?(E F )] the line 
width functions, G r (E F ) = [G a ] t = 1/[E F - H cc - 
S<=i ^i] the retarded and advanced Greens functions, 
and H cc the Hamiltonian of the dashed-box region which 
includes two green-site (light gray) regions and the red- 
site (dark gray) central region [see Fig. 1(a)]. The re- 
tarded self-energy 1^(E F ) due to the coupling to the 
terminal- i can be calculated numerically. 27 To determine 
the longitudinal and Hall resistances mentioned above, 
we applied a bias V across terminal-1 and terminal-4, 
and the currents in the voltage probes (terminals 2, 3, 5, 
and 6) are set to zero. Then from the Landauer-Biittiker 
formula, the voltages V2, V3, V5 and V§ of the voltage 
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FIG. 2: (Color online) The longitudinal resistances Rt and 
Rb vs. En for the different disorder strength W at the energy 
El = —0.1 [panels (a) and (c)] and —0.2 [panels (b) and 
(d)]. The panels (a) and (b) are for R t and the panels (c) 
and (d) are for Rt. The other parameters are the same as in 
Fig. 1(c) and (d). 



probes and the currents I\ and I4 can be calculated. We 
should have I\ = — J4 = In. Finally, the longitudinal 
and Hall resistances are given by R t = (V2 — Vs)/Ii4, 
R b = (Vg - V 6 )/I 1A , Ri = (V 2 - V 6 )/I u and R r = 
(V 3 -V 5 )/I u . 

The recursive Green's function technique 2 is used for 
the computation of the transmission coefficient. The 
Fermi energy E F is set equal to zero as the energy ref- 
erence point. The hopping energy t w 2.75 eV is used 
as the energy unit, which corresponds to 3 x 10 4 K. It 
is reasonable to assume zero temperature condition in 
our calculations because the temperature in the exper- 
iment is only of several Kelvin or sub-Kelvin. Taking 
into account the spin degeneracy, we will use h/2e 2 as 
the resistance unit. The corresponding filling factors vl 
and vr are taken as odd integers ( ±1, ±3, ±5, ...) instead 
of even integers (±2, ±6, ±10, ...). The effect of the con- 
stant magnetic field is addressed by appropriate Peierls 
phase 28 : 20 = (3Vi/2)a 2 B/(j) , where a w 0.142 nm is 
the distance between two neighboring carbons. We will 
take <f> = 0.01 w hich corresponds to the magnetic length 
Ib = y/ft/ieB) « 1.6 nm. 29 We will consider the case 
N = 50 and M = 20, where the area of the central 
region is 21.2 x 5.2 nm 2 and the width of the voltage ter- 
minals is 12.3 nm. The reason for using a width of the 
graphene ribbon far larger than the magnetic length is 
that edge states can not mix except near the boundary 
of the junction. Finally, the disorder is averaged over 
2000 random configurations except in Fig. 4(a) and (b), 
where 500 configurations are taken. 
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III. Numerical results and analysis 

We first study the resistances in the clean graphene 
junction under a strong magnetic field with <j> = 0.01. In 
Fig. 1 (c) and (d), the Hall and longitudinal resistances 
(Ri, R r , R t , and R b ) versus the Dirac energy Er of the 
right region are shown. The Dirac energy El of the left 
region is fixed at -0.1 [Fig. 1(c)] or -0.2 [Fig. 1(d)], 
which corresponds to vi, = 1 and 3, respectively. As 
usual, the Hall resistances Ri and R r display quantized 
plateaus with plateau values ±1, ± 1/3, ± 1/5, ... [in 
units of h/2e 2 ), and Hall plateau of a region only de- 
pends on the filling factor of the region, Ri = 1/vl and 
R r = 1/vr. Furthermore, the Hall resistance plateaus 
are found to be unaffected by the sizes of the central 
region and the presence of the disorder in the central re- 
gion, because that the Hall effect is very robust. We will 
then focus our study on the longitudinal resistances Rt 
and Rb- 

In the presence of the magnetic field, the longitudi- 
nal resistance R t of the top edge usually is not equal to 
Rb of the bottom edge regardless of the junction type 
and disorder strength W, because Hall edge states have 
the chirality which breaks the symmetry of the top and 
bottom edges. With Er < 0, both R t and Rb of the 
clean n-n junction device exhibit perfect plateau struc- 
tures [Fig. 1(c) and (d)]. By considering the carrier trans- 
port along the edges, the plateau values can be analyti- 
cally obtained: 
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for the n-n + regime (0 < vl < vr) and 



R t = and Rb 



1 1 

\vr\ Wl\J 2e 



(3) 



(4) 



for the n + -n regime (0 < ur < vl). The plateau values 
can be understood from the following simple argument. 
It is well known that, for the n-n regime the carriers 
are electron-like and they move clockwise. For the case 
< i>l < vr (the n-n + regime), all edge states at the 
bottom edge are from the right reservoir (i.e. terminal- 
4), so the voltages V§,V§ and V4 are all equal and this 
leads to R b = (V 6 -V 5 )/I u = 0. According to Eq. (2), R t 
is then equal to (1/|^l| — 1 / \v R\)h / 2e 2 . The numerical 
results in Fig. 1(c) and (d) are consistent with the plateau 
values given in Eqs.(3) and (4). For example, in Fig. 1(c) 
where vl = 1, the resistance Rb is zero and Rt can be 
equal to 0, (2/3)h/2e 2 , (4/5)h/2e 2 , ... for v R = 1, 3, 5, 
.... And in Fig. 1(d) where v L = 3, R b is or (2/i)h/2e 2 
and R t can be 0, (2/15)h/2e 2 , (4/21)h/2e 2 , and so on. 

When Er > the device becomes a n-p junction. 
There is no plateau structure for R t and Rb, according 
to Fig. 1(c) and (d). Both R t and Rb are very large and 
they are almost equal. Furthermore, the smaller the fill- 
ing factors are, the larger the longitudinal resistances are. 
For some values of Er, R t and Rb can be over 1000/i/2e 2 
(about 13MS1). This is because for the case of clean n-p 
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FIG. 3: (Color online) The 2-D plots of the longitudinal 
resistance R t (a) and Rt (b) as a function of V g and V c , with 
the parameters N = 50, M = 20, W = 2 and <j> = 0.01. (c) 
The one-dimensional slices shown in (a) (gray solid) and (b) 
(black solid) are plotted to show the dependence of Rt and 
R b on V g with V c fixed at 0.015 and W = 1. 



junction, the edge states in the left and right regions have 
different chiralities and they are well separated in space. 
Hence, edge states mixture can not occur and this leads 
to very large longitudinal resistances. 

Next, we shall discuss how R t and Rb are affected by 
disorder. Fig. 2 shows the dependence of R t and Rb on 
the Dirac energy Er for different disorder strength W. 
E L is fixed at -0.1 (v L = 1) or -0.2 (v L = 3). For 
n-n regime (Er < 0), the plateau structures of Rt and 
Rb still exist for weak and moderate disorder strength. 
And the plateau values are the same as that of the clean 
graphene junction. However, in the case of strong disor- 
der (e.g. W > 3), R t and Rb increase and the plateau 
structures disappear. This is expected as Hall edge states 
begin to be destroyed by strong disorder. For n-p regime 
(Er > 0), R t and Rb are strongly reduced even when 
the disorder is weak. For example, when W = 0.1 (very 
weak) R t and Rb are smaller than 3h/2e 2 for any Er 
(Fig. 2), though R t and Rb can be over 1000/i/2e 2 at 
some values of Er for clean junction [Fig. 1(c)]. This sig- 
nificant decrease results from that the electron-like and 
hole-like Hall edge states start to mix in the vicinity of 
the n-p interface. It should be obvious from our numeri- 
cal result that the top-edge and bottom-edge resistances 
Rt and Rb are not equal and Rb can be negative at some 
specific values of the parameters [Fig. 2(c)]. For suit- 
able disorder strength, the full edge-state mixture occurs 
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FIG. 4: (Color online) The disorder dependence of longitu- 
dinal resistances R t and R), with selected (vl,vr) [which are 
represented by the solid dots in Fig. 3(b)]. (a) and (b) arc 
for Rt and (c) and (d) are for Rb- The other parameters are 
the same as in Fig. 1(c). 



and hence Rt and Rb exhibit plateau structures (see the 
curves of W = 2 in Fig. 2). According to the Landauer- 
Biittiker formula under the condition of full edge-state 
mixture in the central region (junction), the plateau val- 
ues can be analytically obtained: 



Rt = (7^7 + 7^-r)A and R b = 0. (5) 
\vr\ 2e z 



1 



In Fig. 2, the plateau values for some low filling factors 
{ v l,vr) have been labeled and they are well consistent 
with the numerical results. 

In the following we shall base on our model to simulate 
the recent experimental results. In Ref. 17, a difference 
of carrier concentrations between the left and right re- 
gions was introduced by chemical doping and a global 
gate voltage V g was used to control the carrier concen- 
trations of the whole region. So in our theoretical work, 
we instead introduce two voltages V g and V c to control 
the carrier concentrations. Let n L /R be the carrier con- 
centration in the left /right region, ul and ur are re- 
lated to V g and V c in the following way: n/, cx — V g — V c 
and riR cx —V g + V c . Because of the linear dispersion 
relation of graphene, the carrier concentration n L /R is 
approximatively proportional to sgn(E L /R)E L / R . 30 ' 31 So 
we have 

sgn(E L )El = 0(-V B - V c ),sga(E R )E 2 R = p(-V g + V c ), 

(6) 

where (3 is a constant. We assign j3 to be 1 in order to 
simplify the expressions. 

Figs. 3(a) and (b) are the 2-D plots of R t and Rb as 
functions of V g and V c with the disorder strength W = 2. 
The axes El (or v£) and Er (or i/ R ) shown in Figs. 3(a) 
and 3(b) are determined by Eq.(6). The dotted lines in 



Fig. 3(b) are identical to the curves plotted in Fig. 2. 
From the figure we see that Rt and Rb have the following 
symmetry properties: (i) a mirror symmetry with respect 



to the line V g = 



Rb{vL, Vr) 



Rb{—VRi 



-vl) 
-vl), 



(7) 

(8) 



which reflects the inversion of the edge-state chiralities; 
(ii) an inversion symmetry 



Rt{vL,VR) = Rb(vR, v£) 



(9) 



because of the interchange of Rt and Rb by rotating the 
angle 7r round the center of the device. Furthermore, 
both resistances R t and Rb exhibit the plateaus in the 
whole space of the parameters V g and V c . Rt and Rb are 
approximatively constants at fixed filling factors {vl, Vr)- 
But as (i>l, vr) varies, the jump of Rt and Rb occurs so 
that the borders of the filling factors {vl, vr) are clearly 
seen in Fig. 3(a) and (b). Here the resistance plateau val- 
ues in n-n and n-p regimes are coincidental with Eqs.(3), 
(4) and (5). 

In Fig. 3(c), the variations of Rt and Rb along the hor- 
izontal solid lines in Fig. 3(a) and (b) are shown. With 
the range of the gate voltage V g from —0.1 to 0.1, the cor- 
responding filling factors {ul, vr) are (—3, —5), (—1, —3), 
(1, —1), (3, 1), and (5, 3). The figure shows that in unipo- 
lar (n-n and p-p) regime, the top-edge resistance R t is 
always equal to zero. But the bottom-edge resistance 
R b is {2/3)h/2e 2 at (v L ,v R ) = (-1,-3) and (3,1), and 
(2/15)/i/2e 2 at {v l ,vr) = (-3,-5) and (5,3). For the 
bipolar (n-p) regime, {vl,vr) = (1,-1), the plateau of 
R t is 2h/2e 2 and Rb is zero. These resistance plateau 
values are well in agreement with the recent experiment 
results [see Fig. 5(c) in Ref. 17]. 

The disorder strength dependence of the longitudinal 
resistances at the selected filling factors {vl,vr) from 
Fig. 3(b) (represented by solid dots) is shown in Fig. 4. 
Because of the relation between R t and Rb given in 
Eq. (2), R t and Rb have similar characteristics with re- 
spect to the disorder strength W even though R t ^ Rb- 
For n-n regime (see the left panels of Fig. 4), R t and Rb 
in the clean junction (W — 0) have plateau structure. 
The plateau structure still exists when W increases from 
zero. However, when W increases beyond a critical value 
W c , Rt and Rb start to increase from the plateau values. 
The critical value W c is found to be dependent on the fill- 
ing factors {v l ,vr)\ W c w 2.0 for {v L -,vr) = (1,1), (1,3) 
and (1,5), and W c ss 1.5 for (i>l,vr) = (3,3). However, 
for n-p regime (see the right panels of Fig. 4) , R t and Rb 
are very large in the clean junction. When W increase 
from zero, Rt and Rb decrease sharply. For example, 
when W increases from zero to 0.1, Rt is reduced by two 
or three orders of magnitude to a finite value, while Rb 
is reduced to zero. As W continues to increase, both R t 
and Rb decrease to certain plateau values. The plateau 
exists in a certain range of disorder strength, where the 
full mixing of the electron-like and hole-like edge state 
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FIG. 5: (Color online) The resistances R± (a) and Ri, (b) 
vs. the disorder strength W for the different length M of the 
center region. The parameters are El = —0.1, Er — 0.2, 
A = 0.01, and N = 50. 



occurs. For the lowest filling factors {vl,vr) = (1,-1), 
the plateau values Rt = h/e 2 and Rb = exist in a very 
broad range of disorder strength, from 0.1 to 3.0. For 
{vl,Vr) — (1; —3) and (1, —5), the disorder range for the 
existence of resistance plateau are from 1.7 to 2.7 and 
from 1.9 to 2.7, respectively, which is narrower than that 
of (1,-1). For higher filling factors [(vl,vr) — (3,-3) 
or higher], the plateau does not exist, due to the fact 
that it is more difficult to completely mix all edge states 
in the case of high filling factors. This means that the 
resistance plateaus with lower filling factors are easier to 
be observed in experiment. Finally, when W further in- 
crease to be large than 3, Rt and Rb start to increase for 
all cases of filling factors, which indicates that the edge 
states are destroyed by strong disorder. 

Now we study the size effect of the central region on R t 
and Rb- For unipolar graphene junction, the resistances 
Rt and Rb are almost unaffected by the length M and 
width N of central region, except for the case of very 
small N. In the following we focus on the bipolar junc- 
tion case. Fig. 5 shows the dependence of Rt and Rb on 
disorder strength W at different length M for the filling 
factors (vl,Vr) = (1,-3). For other filling factors, the 
results are similar. From Fig. 5, we can see that when 
the central region is short (e.g. M = 10), the full-mixing 
ideal plateaus [Rt — (4/3)/i/2e 2 and Rb = 0] exist in a 
wide disorder range 1.6 < W < 4. This can be under- 
stood as the small size of the central region makes the 
edge states close to each other and hence the full mixing 
of the states occurs. With increasing length M, the full 
mixing is more difficult. Under this condition Rt and Rb 
are enhanced and the corresponding disorder range for 
the ideal plateaus narrows down or even disappears. For 
example, for M — 20, the ideal plateau exists only in 
the disorder range 1.7 < W < 2.7 which is much nar- 
rower than that of M = 10. For larger M such as 50 and 
100, the ideal plateaus do not exist in any case of disor- 
der strength. Although full-mixing plateaus disappear in 
the case of longer junction, extra plateaus [Rt — 2h/2e 2 
and Rb = {2/i)h/2e 2 ] emerge in a long enough junction 



(e.g. M — 100). The existence of the extra plateaus 
can be explained by the partial mixing of edge states. 
The assumption of partial mixing of edge states in long 
enough junctions is reasonable, because part of the Hall 
edge states are near the junction boundary but the other 
Hall edge states are far away from the boundary. Let us 
assume that there are (xl,xr) edge states taking part 
in the full mixing [i.e., the residual — %L, \vr\ — x r) 
edge states in the left and right regions do not participate 
in any mixing] , the plateau values can be analytically ob- 
tained by Landauer-Biittiker formalism: 



Rt 

Rb 



1 



1 



X L Xr 2e z 



1 

xr 



1 



1 h 



(10) 



If xl = \vl \ and xr = \vr\ (i.e. full mixing), the result of 
Eq.(10) reduces to Eq.(5). By taking xl — 1, xr = 1 and 
(vl,vr) — (1,-3), the plateau values given by Eq.(10) 
are R t = 2h/2e 2 and R b = (2/3)h/2e 2 , which are well 
consistent with the numerical results in Fig. 5. 

Finally, we simply discuss how the width N of a 
graphene junction affects Rt and Rb- With increasing 
N at fixed length M, the resistances R t and Rb decrease, 
the resistance plateaus exist in broader disorder range re- 
gardless of the type of the junction (unipolar or bipolar). 
Also, for the large N, some high-filling-factors resistance 
plateaus emerge. This is expected as Hall edge states in a 
wider junction have more chance to mix with each other. 

IV. Conclusion 

In summary, we have investigated the longitudinal re- 
sistances at the two edges of a six-terminal graphene 
junction in the presence of a perpendicular magnetic 
field. By considering the presence of disorder in the vicin- 
ity of the junction interface, the longitudinal resistances 
at opposite edges exhibit different plateaus structures. 
It is found that the plateau values are only determined 
by the Landau filling factors. The numerical results are 
in excellent agreement with the recent experiment by 
Lohmann et ai. Furthermore, for unipolar junctions, 
resistance plateaus exist in clean junctions. The plateau 
structure can be kept in the presence of weak and moder- 
ate disorder, and they are destroyed by very strong dis- 
order. However, for bipolar junctions, the longitudinal 
resistances are very large and do not have any plateau 
structure in the clean case. In the presence of disorder, 
the resistances sharply drop even in the case of very weak 
disorder and they exhibit plateau structures for suitable 
disorder strength. In addition, we study the effect of 
the size of a graphene junction on the resistances and 
find that some extra resistance plateaus emerge in long 
graphene junctions. We explain this by proposing that 
only part of edge states participate in the mixing. 
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